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obtained  by  dividing  computed  potential  jump  by  body  length, 
assuming  that  the  body  is  slender  or  flat  in  the  direction  of 
motion.  The  potential  jump  is  expressed  explicitly  in  terms  of 
the  effective  volume,  i.e.,  the  sum  of  the  displaced  volume  and 
added  mass/density  of  the  submerged  body,  and  the  depth  Froude 
number,  if  a free  surface  is  present,  i^s  a test  of  the  present  speed 
correction  formula,  two  cases  are  considered:  (1)  the  Wigley  para- 
bolic ship  model,  tested  in  both  a small  and  a large  towing  tank, 

(2)  a body  of  revolution  (prolate  spheroid)  tested  in  a 
circular  wind  tunnel.  In  each  case  the  mean-speed  increment 
averaged  over  the  entire  body  surface  is  computed  by  a three- 
dimensional,  finite-element  method  applicable  to  free-surface 
flow  problems.  These  are  shown  to  be  in  good  agreement  with 
those  obtained  by  the  approximate  speed  correction  formula. 

At  high  values  of  Froude  numbers,  the  main  difference  in  the 
total  resistance  coefficients  measured  in  the  two  towing  tanks 
by  Tamura  is  due  primarily  to  difference  in  model  wave  resistance 
computed  for  the  two  tanks  by  a full-fledged,  three-dimensional, 
finite-element  method.  Results  are  also  compared  to  those  ob- 
tained by  using  the  speed  correction  formula  of  Lock  and  Johansen. 

The  present  formula  renders  a better  approximation  than  that  of 
Lock  and  Johansen  when  the  cross  sectional  area  of  a flow  tunnel 
is  not  much  larger  than  the  maximum  cross  section  area  of  the 
body. 
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NOTATION 


A Cross  sectional  area  of  towing  tank 

a One-half  of  length  of  prolate  spheroid 

B Beam  of  ship 

b Radius  of  maximum  cross  section  of  spheroid 

C_  Block  coefficient 

D 

C„  Frictional  resistance  coefficient 

F 

£ Frictional  Resistance  Coefficient 

F 

C^  Total  resistance  coefficient 

C,,  Wave  resistance  coefficient 

W 

C Prismatic  coefficient 

P 

F Hull  fineness  parameter 

F Water  depth  Froude  number 

H 

F Ship  length  Froude  number 

la 

g Acceleration  of  gravity 

H Water  depth 

K Potential  jump  due  to  blockage 

k Partial  form  factor 

Li 

L Body  length;  length  between  perpendiculars 

L Length  of  waterline 

w 

m'  Added  mass  of  submerged  body 

m Added  mass  coefficient 
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«w 

R 

n 

R 

o 
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T 
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Radius  in  cylindrical  coordinates 
Frictional  resistance 

Total  resistance 

Wave  resistance 
Reynolds  number 
Radius  of  tunnel  wall 
Wetted  surface 
Draft  of  ship 

Uniform  incoming  stream  velocity  at  upstream 
infinity 

Mean  speed  due  to  blockage 


u 

o 

W 

x,y,z 

Au 

AG 

v 

P 


Mean  speed  on  the  body  in  unbounded  water 

Width  of  towing  tank 

Right-handed  rectangular  coordinates 

Speed  increment 

1 

Mean-speed  increment  averaged  over  body 
Kinematic  viscosity  of  water 
Density  of  water 


T = (i1,T2,T3)  Tangential  unit  vector 

$ Total  velocity  potential 

$ Total  velocity  potential  in  absence  of 

tank  (or  tunnel)  walls 

4>  Perturbation  velocity  potential 

<J>  Perturbation  velocity  potential  in  absence  of 

tank  (or  tunnel)  wall 
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¥ Displaced  volume;  volume 

V Gradient  operator 
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ABSTRACT 


It  has  recently  been  shown  that  a jump  in  velocity  potential 
exists  between  infinite  upstream  and  downstream  directions  when  a 
body  translates  unformly  along  a channel  of  finite  cross  section 
such  as  a towing  tank  or  wind  tunnel.  In  this  report  a new  block- 
age correction  formula  for  body  speed  is  proposed.  The  speed 
correction  formula  due  to  blockage  is  obtained  by  dividing  com- 
puted potential  jump  by  body  length,  assuming  that  the  body  is 
slender  or  flat  in  the  direction  of  motion.  The  potential  jump 
is  expressed  explicitly  in  terms  of  the  effective  volume,  i.e., 
the  sum  of  the  displaced  volume  and  added  mass/density  of  the 
submerged  body,  and  the  depth  Froude  number,  if  a free  surface 
is  present.  As  a test  of  the  present  speed  correction  formula, 
two  cases  are  considered:  (1)  the  Wigley  parabolic  ship  model, 
tested  in  both  a small  and  a large  towing  tank,  (2)  a body  of 
revolution  (prolate  spheroid)  tested  in  a circular  wind  tunnel. 

In  each  case  the  mean-speed  increment  averaged  over  the  entire 
body  surface  is  computed  by  a three-dimensional,  finite-element 
method  applicable  to  free-surface  flow  problems.  These  are  shown 
to  be  in  good  agreement  with  those  obtained  by  the  approximate 
speed  correction  formula.  At  high  values  of  Froude  numbers,  the 
main  difference  in  the  total  resistance  coefficients  measured  in 
the  two  towing  tanks  by  Tamura  is  due  primarily  to  difference  in 
model  wave  resistance  computed  for  the  two  tanks  by  a full- 
fledged,  three-dimensional,  finite-element  method.  Results  are 
also  compared  to  those  obtained  by  using  the  speed  correction 
formula  of  Lock  and  Johansen.  The  present  formula  renders  a 
better  approximation  than  that  of  Lock  and  Johansen  when  the 
cross  sectional  area  of  a flow  tunnel  is  not  much  larger  than 
the  maximum  cross  section  area  of  the  body. 


ADMINISTRATIVE  INFORMATION 

This  work  was  authorized  and  funded  by  the  Independent  Research  Pro- 
gram at  the  David  W.  Taylor  Naval  Ship  Research  and  Development  Center, 

Task  Area  ZR  011  0201,  Element  61152N. 

INTRODUCTION 

Many  authors  have  investigated  blockage  effect  and  proposed  approxi- 
mate blockage  formulas  to  account  for  towing  tank  or  wind  tunnel  bound- 
1-9* 

aries.  The  first  approximation  concerning  towing  tank  blockage  effects 

date  back  more  than  four  decades.  Owing  to  the  difficulty  encountered  in 
computing  flow  separation,  wake  flow,  free-surface  effects,  etc.,  the  exact 


*A  complete  listing  of  references  is  given  on  page  21. 
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magnitude  of  the  blockage  effect  on  fluid  force  acting  on  a body  is  Coo 
complicated  to  analyze  by  purely  theoretical  means.  However,  these  diffi- 
culties did  not  stop  engineers  from  attempting  to  make  simple  engineering 
approximations  of  the  blockade  problem.  For  engineering  purposes,  compu- 
tation of  a mean-speed  increment  on  a body  due  to  blockage  effects  has  been 
the  main  focus  of  interest  in  order  to  make  a blockage  correction  to 
frictional  drag.  In  the  computation,  the  incremented  change  in  frictional 
drag  due  to  blockage  is  determined  directly  from  the  computed  incremental 
increase  of  mean  speed  over  the  body  surface  caused  by  flow  blockage. 

Two  basic  inviscid  flow-theory  approaches  have  been  previously  em- 
ployed. The  first  approach  is  based  on  the  so-called  one-dimensional, 

mean-flow  theory,  using  the  Kreitner  equation,  which  was  first  obtained  by 
3 

Kreitner  from  the  Bernoulli  and  the  mass  continuity  equations  under  the 

assumption  that  velocity  is  uniform  in  each  cross  sectional  plane.  To 
4 5 

name  a few,  Hughes  and  Kim  used  this  approach.  The  second  approach  is 
based  on  successive  reflection  of  images  in  the  walls  of  a rectangular  tank 
or  simpler  axisymmetric  singularities  in  case  of  axisymmetric  flows.  In 
this  approach,  the  velocity  potential  of  the  flow  inside  a specified  tank 
boundary  can  be  computed  exactly  in  principle;  usually,  the  potential  is 

represented  by  a series  expansion,  and  only  the  first  few  terms  are 

6 12  7 

computed.  Ogiwara,  Tamura,  ’ and  Landweber  and  Nakayama  have  used  the 

latter  approach. 

In  all,  there  exist  about  a dozen  formulas  proposed  for  blockage 

corrections,  and  each  is  somewhat  different  from  the  other.  Some  formulas 

introduce  empirical  correction  factors,  whereas  others  claim  to  be  based 

on  analytical  derivations.  Some  formulas  are  proposed  to  be  used  only  for 

frictional  resistance  corrections,  whereas  other  formulas  are  used  for 

total  resistance  corrections.  An  extensive  review  of  the  subject  has  been 

9 

made  by  Gross  and  Watanabe. 

In  the  present  preliminary  study,  skepticism  is  exercised  about 
proposals  in  speed  correction  formulas  that  can  be  used  to  correct  the 
total  resistance  which  include  the  wave  resistance  in  water  of  finite  depth 


2 


i 


with  sidewalls.  Herein  is  proposed  a new  speed  correction  formula  to  be 
used  only  for  frictional  resistance.  The  wave  resistance  which  has  been 
computed  for  any  towing  tank  and/or  model  conditions  by  using  the  localized 
finite-element  method  previously  developed  by  the  author. ^ ^ It  seems 
to  be  impossible  to  make  blockage  corrections  to  total  resistance  by  using 
only  a single-speed  correction  formula,  even  though  such  formulas  have 
been  proposed  in  the  past. 

The  approach  used  to  derive  a blockage  correction  formula  herein  is 
different  from  the  two  inviscid  flow-theory  approaches  described  pre- 
viously. Derivation  of  a mean-speed  correction  formula  in  this  report  is 
based  on  the  potential  jump  occurring  in  the  three-dimensional  flow  in  a 
towing  tank  or  wind  tunnel.  To  test  the  new  speed  correction  formula, 
numerical  computations  for  a full-fledged,  three-dimensional  wave  resist- 
ance problem  were  made.  The  numerical  mean-speed  increment  on  a spheroid 
was  computed  exactly  for  a circular  wind  tunnel  and  compared  wi th  she 
results  obtained  by  the  new  formula;  results  obtained  agree  reasonably 
well  with  exact  numerical  results. 


BLOCKAGE  CORRECTION 


EXACT  MEAN-SPEED  INCREMENT 


Steady  uniform  flow  past  a ship  fixed  in  a channel  has  been  con- 
12 

sidered;  see  Bai.  The  coordinate  system  is  right  handed  and  rectangular. 
Under  the  usual  assumptions,  steady  uniform  flow  may  be  described  by  a 
total  velocity  potential  0 defined  by 

<P  (x , y , z)  = Ux  + <j>  (x,  y,  z)  (1) 

where  is  the  perturbation-velocity  potential  in  a channel  of  finite 
cross  section.  Similarly  the  total  velocity  potential 


$o  (x , y , z)  = Ux  + </>o  (x,  y,  z) 


3 


r 

I 


is  defined  to  describe  the  flow  about  the  same  body  in  an  unbounded  fluid, 
i.e.,  in  the  absence  of  channel  boundaries.  The  fluid  speed  on  a body 
surface  in  general  increases  due  to  the  blockage  effect  when  compared  with 
that  of  unbounded  fluid.  However,  the  speed  increment  on  the  body  surface 
is  not  uniform  over  the  entire  surface.  For  example,  the  forward  stag- 
nation point  of  an  axisymmetric  body  remains  the  same  whether  in  an  un- 
bounded fluid  or  in  a wind  tunnel  of  circular  cross  section.  Nevertheless, 
a mean  speed  correction  has  been  traditionally  employed  for  the  blockage 
correction  mainly  due  to  its  simplicity.  To  describe  a mean-speed  incre- 
ment, speed  increment  due  to  blockage  locally  on  the  body  surface  is 
defined  as 


Au  = V ($-$  ) • x 
o 

= V (Q-Q  ) • X 
o 


(3) 


where  x = (X^, 12,13)  is  a unit  tangential  vector  on  the  body  surface;  x^ 
is  the  component  along  the  x-axis,  i.e.,  the  longitudinal  direction,  and 
X?  and  x^  are,  respectively,  the  normal  and  tangential  components  in  the 
cross  sectional  plane  of  the  body.  Then  the  "exact"  mean  speed  increment 
averaged  over  the  entire  submerged  body  surface  is  given  by 


Au 


1 

S 


V ((jj-cj)^)  • x ds 


(4) 


where  S is  the  wetted  surface  area,  and  x is  specified.  One  natural  way 
of  specifying  x would  be  as  the  unit  potential  flow  streamline  vector  on 
the  body.  However,  streamlines  on  a body  in  bounded  and  unbounded  flows, 
described  by  $ and  respectively,  do  not  coincide  in  general,  except  in 
the  special  case  of  an  axisymmetric  body  in  a flow  facility  of  circular 
cross  section.  In  the  case  of  a ship  hull,  if  x = (1,0,0),  and  = 

2 • L • T under  the  assumption  that  the  ship  is  thin,  Equation  (4)  can  be 
reduced  to 


1 


4 


(5) 


Au  = u - 


u 

o 


where 


0 

“ = irKj  [*  (2  *y*°)  - * (-  \ >y*°)  ]dy 

0 

% = rH  J [♦o  (1  >y>°)  " t (~  2 *y*°)  ] dy 

-T 

where  L and  T are  the  ship  length  and  draft,  respectively.  In  Equation  (5), 
the  draft  T is  assumed  to  be  uniform  from  the  bow  at  x =-  L/2  to  the  stern 
at  x = L/2;  the  centerplane  of  the  ship  is  on  z = 0. 

Similarly,  for  a slender  axisymmetric  body  of  revolution  in  a wind 
tunnel  of  circular  cross  section,  the  mean-speed  increment  averaged  over 
the  body  surface  is  given  by 


Au  = u - u 

o 


where 


u 


= \ , R=0 
=-  \ , R=0 


u 

o 


X = | , R=0 
x =-  ~ , R=0 


(6) 


where 

and  the  peripheral  length  along  a body  meridian  is  approximated  by  the 
body  length,  assuming  that  the  body  is  slender. 
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APPROXIMATE  MEAN-SPEED  INCREMENT 


In  this  subsection,  the  method  of  obtaining  an  approximate  speed 

correction  formula  is  given,  based  on  the  potential  jump  discussed  earlier 
12 

by  Bai.  Define  K as  a jump  in  the  velocity  potential  <J>  given  in 
Equation  (1)  between  the  infinite  upstream  and  downstream  directions.  The 
potential  jump  K is  given  by  integrating  the  speed  increment  along  a line 
in  the  fluid  from  a point  infinitely  far  upstream  to  a point  infinitely 
far  downstream.  Numerical  solutions  for  practical  ship  forms  at  sub- 
critical  speeds  in  towing  tanks  and  for  slender  bodies  of  revolution  in 
wind  tunnels  indicate  that  most  of  the  potential  jump  occurs  along  the 
body  length.  This  finding,  observed  in  numerical  solutions,  will  be  used 
as  the  basis  for  obtaining  the  present  approximate  formula  for  the  speed 
correction.  It  is  possible  to  prove  this  empirical  finding  by  showing 
that  the  values  of  the  potential  at  the  upstream  and  downstream  stagnation 
points  are  approximately  equal  to  the  corresponding  asymptotic  values  of 
the  potential  in  the  simple  case  of  axisymmetric  flow.  However,  the  proof 
will  not  be  discussed  here.  Thus,  the  mean  speed  increment  u due  to 
blockage  is  approximated  by 


4J  - £ (7) 

In  a recent  simple  analysis,*  the  expressions  for  the  potential  jump 

K in  terms  of  the  effective  volume  and  the  depth  Froude  number  F in  three 

H 

dimensions  with  a free  surface  were 


= (V  -Hn'/p)  U 
WH  (1-Fh2) 


*A  more  detailed  analysis  in  general  cases  has  been  submitted  in  a 
paper  to  the  Journal  of  Fluid  Mechanics  (1978). 
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where  ¥ = displaced  volume 

p = density  of  water 
U = towing  speed 
W = tank  depth 
H = water  depth 

m'  = m'(Fjj)  = added  mass  in  the  longitudinal  direction 
Fh  = U/ /gH 

In  the  derivation  of  Equation  (8),  it  is  assumed  that  the  waterplane  area 
of  the  ship  hull  is  so  thin  that  a line  integral  term  is  neglected,  and  the 
body  boundary  condition  is  satisfied  exactly  on  the  body  surface.  From 
Equations  (7)  and  (8)  is  obtained 


Au 

U 


¥ + m1  /p 
AL  U-Fh2) 


(9) 


where  A = WH  is  the  cross  sectional  area  of  the  tank.  It  is  of  interest 
to  note  that  when  the  value  of  g approaches  infinity,  F approaches  zero, 

ri 

and  Equation  (9)  reduces  to  the  case  of  a wind  tunnel,  where 


Au  _ ¥ + m'  /p 
U " AL 


(10) 


It  is  also  of  interest  to  note  that  when  the  body  boundary  condition  is 
linearized,  i.e.,  satisfied  on  the  body  centerplane,  Equations  (9)  and  (10) 
further  reduce  to 


Au 

U 


¥ 

AL  (1~Fh2) 


(9’) 


in  the  presence  of  a free  surface  and 
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(10') 


Au  _ V_ 

U = AL 

in  the  absence  of  a free  surface. 

APPLICATIONS 

TOWING  TANK  EXPERIMENT 

To  test  the  new  blockage  correction  formula,  three  sets  of  compu- 
tations were  first  made  for  the  same  model  in  three  different  towing  tanks. 

1 2 

The  first  two  tanks  had  the  dimensions  given  by  Tamura;  ’ see  Table  1. 

The  third  tank  was  approximately  four  times  greater  in  cross  sectional  area 
than  the  large  tank  listed  in  Table  1,  i.e.,  W = 24  m and  H = 12  m.  The 
specific  ship  model  considered  was  the  Wigley  parabolic  model  (Model  M1719 
in  Tamura),  and  the  equation  of  the  hull  surface  was  given  by 


where  L/B  = 10,  and  T/L  = 0.0625.  The  geometric  particulars  of  the  models 
have  been  given  in  Table  2. 

In  the  computations,  the  ship  hull  boundary  condition  was  linearized; 

thus,  speed  correction  formula  (Equation  (9'))  was  used.  To  test  the 

present  mean-speed  correction  formula,  computations  were  also  made  from 

Equation  (5)  the  exact  mean-speed  increment  averaged  over  the  hull  surface 

. 12 

from  the  local  velocities  obtained  by  the  finite-element  method.  In 
computing  the  value  of  uq  from  Equation  (5),  the  numerical  result  for  the 
extra  large  tank  was  used  in  place  of  the  perturbation  potential  for  un- 
bounded water  d>  because  the  effect  of  the  tank  wall  and  the  bottom  was 
o 

found  to  be  negligibly  small.  Comparisons  between  the  "exact"  and 
approximate  mean-speed  increments  are  given  in  Table  3.  Agreement  is 
reasonably  good.  It  should  be  noted  in  Table  3 that  the  exact  mean  speed 
averaged  on  the  hull  surface  u^,  defined  by  Equation  (5),  is  not  only 
nonzero  but  also  independent  of  Froude  number.  It  should  also  be  noted 
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TABLE  1 - DIMENSIONS  OF  SMALL,  LARCE,  AND  EXTRA  LARGE 
TOWING  TANKS1,2 


Width  in  meters 

6.09 

12.5 

24 

Mean  Water  Depth  in  meters 

3.555 

6.268 

12 

TABLE  2 - WIGLEY  PARABOLIC  MODEL  (TAMURA  MODEL  M1719) 


Length  between  Perpendiculars  in  meters 

8.000 

Length  of  Waterline  in  meters 

7.984 

Beam  in  meters 

0.800 

Draft  in  meters 

0.500 

Volume  in  meters 

1.422 

Wetted  Surface  in  meters 

9.408 

Block  Coefficient 

0.4453 

Prismatic  Coefficient 

0.6680 

TABLE  3 - COMPARISONS  OF  MEAN-SPEED  INCREMENT,  COMPUTED  BY 
NUMERICAL  RESULTS  AND  BY  PRESENT  FORMULA  FOR 
F = 0.4* 


Exact  Numerical  Results 

Equation  (9') 

Tank 

u /U 
o 

u/U 

Au/U  = (u-u  )/U 
o 

Au/U 

Small 

0.017198 

0.030514 

0.0133 

0.0128 

Large 

0.017198 

0.019425 

0.0022 

0.0029 

*Results  of  extra  large  tank  were  used  to  compute  u 
as  discussed  in  text. 


that  the  free  surface  effect  on  the  velocity  profile  on  the  body  surface 
would  be  significantly  dependent  upon  whether  the  hull  is  in  a shallow 
towing  tank  or  in  unbounded  water.  The  present  study  indicates  that  the 
approximate  speed  correction  formula  satisfactorily  treats  the  seemingly 
complicated  free-surface  effect  on  the  mean-speed  increment  on  the  body. 

The  total  resistance  coefficient  C^,  determined  experimentally  by 
Tamura,  and  the  wave  resistance  coefficients  computed  by  the  finite 
element  method,  are  given  in  Table  4. 

In  presenting  our  results,  the  total  resistance  coefficient  C T and 

the  wave  resistance  coefficient  C are  defined  as 

w 


CT  = RT/f  1)2  *2‘ ' 3 

C = R /f  U2  ¥2/3 
w w/  2 


(12) 


where  R^  and  R^  are,  respectively,  the  total  and  wave  resistances. 

A 

The  frictional  resistance  coefficients,  C_,  and  C_,  are  defined  by 

r r 


c 


F 


S 

o 


A 


C 


F 


r/f  u 


2 ¥2/3 


(13) 


where  Sq  is  the  model  wetted  surface  area.  The  model  length  Froude  number 

F and  Reynolds  number  R are  defined  by 
L n 


Fl  = U//gL 

R = V//UL 
n 


(14) 
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where  v is  the  kinematic  viscosity  of  water.  Here  the  Reynolds  number  R^ 

is  obtained  by  assuming  that  the  freshwater  temperature  in  two  towing 

tanks  was  20  C.  Table  A results  are  given  in  Figure  1;  the  wave  resistance 

computed  for  the  extra  large  tank  was  taken  to  be  the  same  as  for  unbounded 

water,  already  mentioned.  In  Figure  1,  hull  wave  resistance  in  the  large 

tank  is  very  close  to  that  for  the  extra  large  tank.  Thus,  the  blockage 

effect  on  wave  resistance  is  very  small  for  the  large  tank.  Also,  the 

main  difference  in  the  total  resistance  coefficients  measured  in  the 

small  and  large  towing  tanks  is  due  primarily  to  the  difference  in  the 

model  wave  resistance  computed  for  the  two  tanks. 

Table  A gives  the  speed  corrections  computed  from  Equation  (9')  along 

with  the  corrected  values  of  (C_-C  ).  The  corrected  value  of  (C  -C  ) is 

2 ” w T w 

given  by  (C„-C  ) (U/(U+Au))  . Table  5 gives  the  frictional  resistance  co- 

^ w A 

efficients  C and  C , computed  from  International  Towing  Tank  Conference 
r r 

(ITTC)  (1957)  and  American  Towing  Tank  Conference  (ATTC)  friction  formulas. 

In  the  present  study,  it  is  assumed  that  the  total  resistance  less  the 

computed  theoretical  wave  resistance  is  approximately  equal  to  the 

frictional  resistance,  since  the  ship  hull  is  thin  and  smooth,  i.e.,  form 

13 

drag  is  assumed  to  be  negligibly  small.  If  we  make  use  of  the  Granville 
correlation  of  partial  form  factor  k^  with  hull-fineness  parameter 
F = C_  / (B/L) (2T/L)  for  the  Wigley  parabolic  model  with  F = 0.5,  we  find 
that  k^  = 0.0A;  i.e.,  form  drag  is  estimated  to  be  only  A percent  of  the 
frictional  drag  and  a still  lower  percentage  of  the  total  drag.  Accord- 
ingly, speed  correction  Equation  (9')  was  applied  to  the  resistance 

component  C = C„  - C to  correct  for  blockage  effect.  Results  given  in 
r r W 

Table  5 are  shown  in  Figure  2.  In  Figure  2,  the  corrected  values  of 
(C„-C  ) are  lower  than  the  values  of  C„  given  by  ITTC  and  ATTC  friction 
formulations,  indicating  negative  hull  form  drag,  which  is  not  acceptable. 
In  other  words,  if  the  form  drag  coefficient  and  other  corrections  had 
been  added  to  the  values  of  ITTC  and  ATTC  friction  coefficients,  this 
discrepancy  would  be  even  larger.  The  discrepancy  seems  to  have  been 
caused  by  computed  values  of  the  wave  resistance  being  too  large. 
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TABLE  4 - RESISTANCE  COEFFICIENTS  OF  THE  WIGLEY  PARABOLIC  MODEL  AT 
TWO  DIFFERENT  TOWING  TANKS  (MODEL  M1719) 


fl 

CT 

Experiment 
by  Tamura 

C 

w 

Numerical 
by  Bai 

VCw 

Au/U 

(CT-Cw> 

Corrected 

(fl> 

Corrected 

0.325 

0.0353 

0.00959 

0.0257 

0.0108 

0.0251 

0.329 

Small 

Tank 

0.350 

0.0348 

0.00981 

0.0250 

0.0113 

0.0244 

0.354 

0.375 

0.0397 

0.01543 

0.0243 

0.0120 

0.0237 

0.380 

0.400 

0.0480 

0.02364 

0.0244 

0.0128 

0.0237 

0.405 

0.425 

0.03106 

0.325 

0.0344 

0.00971 

0.0247 

0.0026 

0.0246 

0.326 

0.350 

0.0331 

0.00785 

0.0253 

0.0027 

0.0251 

0.351 

Large 

Tank 

0.375 

0.0361 

0.01203 

0.0241 

0.0028 

0.0240 

0.376 

0.400 

0.0414 

0.01865 

0.0228 

0.0029 

0.0226 

0.401 

0.425 

0.02487 

Values  in  both  tanks  are 

converted  to 

the  case  without  blockage. 

TABLE  5 - FRICTIONAL  RESISTANCE  COEFFICIENTS  C_  AND  C_, 

r r 

COMPUTED  FROM  ITTC  (1957)  AND  ATTC  FRICTION  FORMULAS, 
AT  A FRESHWATER  TEMPERATURE  OF  20  C 


ITTC 

1957 

ATTC 

fl 

R 

n 

CF 

A 

CF 

CF 

A 

CF 

0.325 

3.861  x 106 

0.003565 

0.0265 

0.003444 

0.0256 

0.350 

4.158  x io6 

0.003516 

0.0261 

0.003400 

0.0253 

0.375 

4.455  x io6 

0.003470 

0.0258 

0.003360 

0.0250 

0.400 

4.752  x io6 

0.003429 

0.0255 

0.003323 

0.0247 
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In  the  numerical  computation  of  wave  resistance  by  the  finite-element 
method,  44  nodes  on  the  ship  hull  surface,  i.e.,  on  the  centerplane,  and 
1496  nodes  for  the  entire  fluid  domain  were  taken.  One  may  expect  more 
refined  results  by  reducing  the  size  of  finite  elements.  To  treat  low 
values  of  Froude  number  accurately,  smaller  and  more  elements  are  necessary. 

WIND  TUNNEL  EXPERIMENT 

As  a second  example,  the  blockage  effect  was  considered  for  a wind 
tunnel  having  a uniform  circular  cross  section  of  radius  R^.  The  specific 
body  geometry  considered  was  a prolate  spheroid  with  its  meridian  profile 
given  by 


(15) 


for  the  special  case  when  a/b  = 4. 

The  potential  flow  for  the  axisymmetric  boundary  configurations  con- 
sidered herein  could  have  been  computed  by  the  conventional  method  of 

integral  equations;  i.e. , the  axial  source  and  doublet  distributions  or 
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the  vortex  sheet  on  the  surface,  etc.,  as  discussed  in  Landweber. 

However,  the  velocity  potential  has  been  computed  by  the  finite-element 

method.  Computations  have  been  made  for  seven  values  of  RQ/b  = 1.25, 

1.5,  2,  3,  4,  5,  and  15  all  for  a/b  = 4.  When  R^/b  = 15  was  computed,  the 

effect  of  the  tunnel  wall  on  the  body  surface  was  negligibly  small  as  if 

the  body  were  moving  in  an  infinite  fluid.  The  value  of  uq/U  defined  in 

Equation  (5),  computed  by  using  the  result  of  R0/b  = 15,  was  0.08185, 

whereas  that  computed  by  using  the  exact  analytic  result  for  the  unbounded 

water,  i.e.,  R^/b  = 00 , given  in  Lamb^  was  0.08156. 

The  computed  velocity  potential  <}>  is  shown  in  Figure  3 for  three 

values  of  R /b  = 1.25,  1.5,  and  15.  To  illuminate  the  assumption  made  to 

o 

obtain  the  present  approximate  mean-speed  correction.  Figure  3 shows 
straight  lines  drawn  from  the  origin  to  the  asymptotic  values  of  K/2  at  the 
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LEGEND: 

ON  BODY 
AT  R/b  = 1.25 

LINEAR  POTENTIAL  VARIATION 
ASSUMED  IN  THE  PRESENT  SPEED 
CORRECTION  FORMULA  (THE 
SLOPE  IS  THE  SPEED  CORRECTION) 

Figure  3 - Velocity  Potential  for  a Spheroid  (a/b=4) 
in  a Wind  Tunnel  with  a Circular  Cross  Section 
of  Radius  R 


p 


the  downstream  stagnation  point  x = L/2.  The  slope  of  each  straight  line 

is  equal  to  the  speed  correction  defined  by  Equation  (7) . Owing  to  the 

skew  symmetry  of  the  potential  with  respect  to  x = 0,  the  result  for  the 

upstream  half-body  can  be  obtained  from  the  downstream  potential  shown  in 

Figure  3.  The  velocity  potential  increases  monotonically  from  a value 

slightly  lower  than  - K/2  at  the  upstream  stagnation  point  to  a value 

slightly  higher  than  K/2  at  the  downstream  stagnation  point  on  the  body 

surface.  However,  the  potentials  at  R = 1.25b  approach  monotonically  the 

asymptotic  values  at  both  ends  for  R /b  = 1.25  and  R /b  = 1.5. 

o o 

In  Table  6 the  approximate  mean  speed  correction  given  by  Equation 
(10)  is  compared  with  the  exact  mean-speed  correction  computed  from 
Equation  (6).  Table  6 also  gives  the  speed  correction  obtained  by  the 

g 

Lock  and  Johansen  formula,  which  is  given  in  Pope  as 


Au 

U 


(16) 


When  R /b  < 3,  our  approximate  results  show  better  agreement  with  the  exact 
o 

numerical  results  than  with  those  of  Lock  and  Johansen. 

In  Figure  4 computed  values  of  the  added  mass  coefficient  and  the 

mean  speed  correction  Au/U  are  shown  as  a function  of  b/R^.  In  Figure  4, 

note  that  for  b/R  > 0.765,  the  contribution  of  the  added  mass  to  the 
o 

speed  correction  in  Equation  (10)  is  more  dominant  than  the  contribution 
of  the  displaced  volume,  i.e.,  m = m'/P^  > 1.  This  finding  indicates  that 
a crude  blockage  correction,  based  on  only  the  local  cross  sectional  area 
of  the  body  using  one-dimensional  theory,  cannot  always  give  a good 
approximation  of  the  mean-speed  correction  when  the  added  mass  coefficient 
is  not  small. 


CONCLUSIONS 

In  the  present  study  a new  mean-speed  formula  for  corrections  caused 
by  blockage  is  proposed.  The  approximate  formula  is  tested  by  comparing 
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m AND  Au/U 


TABLE  6 - COMPARISONS  OF  MEAN-SPEED  INCREMENTS  ON  A 
SPHEROID  IN  A WIND  TUNNEL,  COMPUTED  BY  A NUMERICAL 

METHOD,  APPROXIMATE  FORMULAS  (u  /U=0. 0813557 

o 

OBTAINED  BY  LAMB  WAS  USED) 


R /b 
o 

u/U 

Au/U 

Exact 

Present 
Formula 
(Equation  (6)) 

Lock  and  Johansen 
(Equation  (16)) 

2.25 

0.  98204 

0.90050 

0.93965 

1.22419 

1.5 

0.52285 

0.44129 

0.47716 

0.70844 

2 

0.26702 

0.18546 

0.21559 

0.29888 

3 

0.14442 

0.06287 

0.08505 

0.08856 

4 

0.11088 

0.02932 

0.04624 

0.03736 

5 

0.09748 

0.01593 

0.02920 

0.01913 

15 

0.08185 

0.00030 

0.00320 

0.00071 

Figure  4 - Added  Mass  Coefficient  m and  Speed  Correction 
Au/U  for  a Spheroid  in  a Circular  Wind  Tunnel 
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* 


it  with  an  exact  numerical  mean-speed  correction,  computed  by  the  finite- 
element  method  for  both  a towing  tank  experiment  and  a wind  tunnel  experi- 
ment, The  two  predictions  are  shown  to  be  in  good  agreement  for  both  \ 

facilities.  It  is  shown  that  the  effect  of  added  mass  coefficient  on  the 

] 

speed  correction  of  a body  is  very  significant  as  the  blockage  effect  in- 
creases. It  is  also  found  that  the  main  difference  in  the  total  resistance 
coefficient  measured  in  a large  and  a small  towing  tank  is  due  primarily 
to  the  difference  in  the  model  wave  resistances  computed  for  the  two  tanks. 

Further  investigation  is  necessary  to  take  into  account  other  blockage 
corrections  due  to  viscous  effects  such  as  flow  separation  and  wake  dis- 
placement thickness  effects. 
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DTNSRDC  ISSUES  THREE  TYRES  OF  REFORTS 


1.  DTNSRDC  REPORTS.  A FORMAL  SERIES.  CONTAIN  INFORMATION  OF  PERMANENT  TECH- 
NICAL VALUE.  THEY  CARRY  A CONSECUTIVE  NUMERICAL  IDENTIFICATION  REGARDLESS  OF 
THEIR  CLASSIFICATION  OR  THE  ORIGINATING  DEPARTMENT. 

2.  DEPARTMENTAL  REPORTS.  A SEMIFORMAL  SERIES.  CONTAIN  INFORMATION  OF  A PRELIM 
INARY,  TEMPORARY,  OR  PROPRIETARY  NATURE  OR  OF  LIMITED  INTEREST  OR  SIGNIFICANCE. 
THEY  CARRY  A DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION. 

3.  TECHNICAL  MEMORANDA.  AN  INFORMAL  SERIES.  CONTAIN  TECHNICAL  DOCUMENTATION 
OF  LIMITED  USE  AND  INTEREST.  THEY  ARE  PRIMARILY  WORKING  PAPERS  INTENOEO  FOR  IN- 
TERNAL USE.  THEY  CARRY  AN  IDENTIFYING  NUMBER  WHICH  INDICATES  THEIR  TYPE  AND  THE 

» NUMERICAL  CODE  OF  THE  ORIGINATING  DEPARTMENT.  ANY  DISTRIBUTION  OUTSIDE  DTNSRDC 

MUST  BE  APPROVED  BY  THE  HEAD  OF  THE  ORIGINATING  DEPARTMENT  ON  A CASE-BY-CASE 
BASIS. 


